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I.  INTRODUCTION 


J 


The  position  of  a  source  of  radiation  is  often  estimated  from  bearing  observations 
made  from  two  or  more  known  locations.  If  we  let  (x,  y)  denote  the  coordinates  of  the 
source  and  a.  its  true  bearing  from  a  station  with  site  location  (X.,  Y.)  for  i  =  1,  2,  . . 

.  ...n,  then,  if  precise  measurements  were  possible,  the  geometrical  relationships  link¬ 
ing  the  known  quantities  X.,  Y.,  o\  and  the  unknown  quantities  x  and  y  could  be  used  to 

determine  the  required  coordinates  of  the  source.  The  geometrical  relationships  can  be 
formulated  for  each  station  separately  and  therefore  have  the  form 

=  f(X.,  Y.,  x,  y)  =  f.(x,  y)  (1-1) 

for  i  =  i,  2, . . .  n  . 

In  practice,  the  angular  observations  will  be  subject  to  errors  of  several  kinds.  Since, 
however,  the  treatment  of  systematic  errors  in  the  measurement  equipment  and  errors 
due  to  spurious  sources  of  radiation  usually  consists  of  procedures  for  refining  the  ob¬ 
served  data  before  estimating  the  fix,  we  shall  take  as  our  starting  point  the  possibly- 
modified  data  set: 

W*!  »-2) 

for  i  =  1,  2  . .  .n, 
or,  in  vector  notation 

9=a+yr  (i.2a) 

where  the  Y.  represent  random  measurement  errors  having  the  joint  probability  density 
function: 

p(^  =  p(*rY2, . Yn)  .  (1.3) 

The  problem  of  fix  estimation  is  one  of  constructing  an  algorithm  which,  for  the  appro¬ 
priate  site  geometry  and  for  an  assumed  error  distribution,  provides  an  optimal  estimate 
of  the  source  position  from  the  data  furnished  by  the  observations. 
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II.  SOME  CRITICISMS  OF  CURRENT  ESTIMATION  PROCEDURES 

If  it  is  accepted  that  the  true  bearings  have  fixed,  but  unknown  values,  then,  since  all 
of  the  pertinent  geometrical  information  is  contained  in  the  nonlinear  relations  (1. 1),  the 
problem  is  seen  to  consist  of  estimating  a  known  nonlinear  function  of  the  parameter  a 
of  the  joint  density  function  p(£;a)  the  form  of  which  is  readily  obtained  from  that  of  the 
assumed  error  density  P(^)-  Consequently,  it  is  not  surprising  that  for  many  systems 
in  current  use,  the  estimating  procedures  are  based  on  the  principle  of  maximum  likeli¬ 
hood.  References  (a)  through  (d)  give  procedures  obtained  by  this  approach. 

The  use  of  maximum  likelihood  methods  for  fix  estimation  has  been  subject  to  some 
criticisms.  These  include  the  following: 

(a)  Several  factors  limit  the  numbers  of  bearing  observations.  There  is  the  avail¬ 
ability  of  suitable  station  sites  and  the  cost  of  measurement  equipment.  Also,  if  sur¬ 
veillance  is  exercised  over  an  extended  area  or  if  the  duration  of  the  radiation  is  small 
enough,  only  some  stations  may  have  observations  to  report  and  of  these,  some  may  be 
rejected  as  being  invalid.  For  what  is  more  often  than  not  a  small  sample  problem, 
the  properties  of  consistency  and  asymptotic  efficiency  are  clearly  irrelevant,  [e  ] 

(b)  Although  it  is  fairly  straightforward  to  formulate  the  likelihood  equations,  the 
nonlinear  character  of  the  geometrical  relationships  (1.1)  generally  precludes  a  solution 
in  closed  form  [f,  c].  This  leads  to  solutions  in  specific  cases  which  are  direct  approxi¬ 
mations  or  are  based  on  iterative  techniques  or  which  combine  such  procedures.  Un¬ 
fortunately,  the  derivations  of  these  solutions  seldom  include  adequate  treatments  of  the 
implications  of  the  approximations  or  the  convergence  properties  of  a  proposed  iterative 
procedure. 

(c)  The  method  of  maximum  likelihood  makes  no  provision  for  making  use  of  prior 
knowledge.  This  is  almost  always  present  in  some  degree.  Thus,  not  only  is  valid  in¬ 
formation  not  incorporated,  in  addition  there  is  the  possibility  that  a  correct  estimate  of 
the  most  probable  position  of  the  source  may  be  in  conflict  with  this  prior  knowledge . 

The  example  of  the  estimated  position  being  on  land  for  an  emitter  known  to  be  shipbome 
is  not  necessarily  a  consequence  of  anomalies  in  either  the  data  or  the  fix  algorithm. 


III.  POSITION  FINDING  USING  BAYES*  PROCEDURES 

If,  in  contrast  to  the  assumptions  of  the  previous  section,  the  coordinates  of  the 
radiating  source  are  assumed  to  be  random  variables  for  which  a  realistic  joint  density 
can  be  assigned  from  knowledge  existing  before  the  observations  are  made,  the  problem 
maybe  formulated  within  a  Bayesian  framework.  This  formulation  permits  the  available  prior 
information  to  be  incorporated  in  theestimation  procedure.  For  a  postulated  prior  density 

A  A 

cr(x,  y),  the  required  Bayes*  estimator  is  given  by  the  values  of  x  and  y  which,  for  each 
possible  set  of  bearing  observations  Q,  minimize  the  a  posteriori  risk,  namely: 


where, 


IM  y;  x,  y)  p(x,  y  I  q)  dx  dy 

P(x.  y|fl>  - 

=  p(fl|x,  y)or(x,  y) 


X,  y)or(x,  y)  dx  dy 


(3.1) 


(3.2) 

(3.2a) 


is  the  a  posteriori  density,  expressed  here  in  terms  of  known  quantities,  and 

L(x,  y, ;  x,  y)  is  an  appropriate  cost  function.  Where  this  is  a  quadratic  function  of  the 

estimation  errors  (x  -  x)  and  (y  -  y),  or  alternatively,  certain  conditions  listed  in 
reference  (g)  are  known  to  hold,  the  Bayes*  estimator  is  simply  the  mean  of  this  a 
posteriori  density  function.  It  may  be  obtained  from  quantities  assumed  to  be  known  by 
means  of  the  two  expressions: 


■//'  p(x,  yle)  dx  dy 


ff*  p(5l  x,  y)  cr(x,  y)  dx  dy 
ffp(8\  x,  y)  cr(x,  y)  dx  dy 


(3.3a) 


rr  i  //y  p  (*|x.  y)  o-(x,  y)  dx  dy 

=  //y  P  (x.  y  0)  dxdy  =  - 

JJ  P  (5|  x,  y)  or(x,  y)  dx  dy 


(3.3b) 


These  standard  forms  represent  but  the  starting  point  for  the  derivation  of  a  practical 
algorithm  for  fix  estimation.  Although  they  make  formal  provision  for  the  inclusion  of 
prior  knowledge  this  will,  in  general,  be  of  an  empirical  character.  Consequently,  any 
algorithm  based  on  the  above  Bayesian  formulation  presupposes  that  the  conversion  of  this 
empirical  knowledge  to  an  acceptable  prior  density  has  been  achieved.  Furthermore, 
there  is  no  evidence  that  this  approach  offers  any  advantage  over  that  based  on  maximum 
likelihood  in  providing  solutions  in  closed  form.  To  illustrate  these  points,  let  us  examine 
the  special  case  most  frequently  considered  in  the  material  referenced.  For  this  case, 
the  assumptions  are  as  follows: 
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(a)  The  surveillance  area  is  small  enough  for  the  curvature  of  the  earth's  surface 
to  be  neglected. 

(b)  For  all  stations  the  observation  errors  do  not  depend  on  their  respective  distances 
from  the  source  of  radiation  and  are  normally  (though  not  necessarily  identically)  distrib¬ 
uted  with  mean  zero. 

(c)  The  observation  errors  for  each  station  are  independent. 


The  first  assumption  determines  the  geometrical  relationship  applicable  to  this  case. 
Assuming  angular  measurements  from  the  direction  indicated  by  a  positive  y  coordinate, 
it  is  readily  shown  that,  for  this  case,  the  relation  (1.1)  has  the  form: 

x  -  X. 

ai  =  (x,  y)  =  arctan  —  y1 —  •  (3.4) 

y 

for  i  =  1,  2  . . .  n,  where  n  is  the  number  of  valid  observations.  The  remaining  as¬ 
sumptions,  together  with  equations  (1.3)  and  (1.2a)  yield  the  following: 


-  1/2  -  1  ‘  \2/2oi2 

p(Yi)  =  (2tt)  *e  1  1 


(3.5 


P(0la)  =  n 

i  =1 


v-  1/2  -  1 

(2tt)  tr.  e 


2ar? 


1  <V"i>2 


(3.6) 


:  (2rr)  *n/2  (o1  .. 


.-I 

a  )  e 


-  1/2  £  7  (6i  ' 


i  =  1 


(3.6a) 


P(0|x,  y)  =  (2tt) 


u 

-n/2  .  ,-l  -1/2  y  V  {9i  -fi  <x.y)}‘ 

(a  ...  er )  e  '  a.  1  ’ 


i=l 


(3.7) 


where,  from  (3.4). 

x  -  X. 

fj(x,  y)  =  arctan  -  — yl —  . 


Substitution  in  (3.2a)  gives  the  a  posteriori  density  corresponding  to  the  arbitrary  prior 
density  o(x,  y),  namely: 
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and  substituting  in  (3.3a)  and  (3.3b)  gives  the  expressions  for  the  corresponding  fix 
estimator,  namely: 


These  expressions  serve  to  illustrate  the  general  problem.  First,  a  realistic  prior 
density  has  to  be  assigned  on  the  basis  of  existing  prior  knowledge,  then,  means  must  be 


found  for  reducing  the  resulting  complex  mathematical  forms  to  feasible  computational 
procedures.  In  the  following  section  a  practical  solution  to  these  problems  will  be 
advanced.  Essentially  it  consists  of  formulating  the  prior  knowledge  in  a  way  which 
simplifies  the  mathematical  relationships  which  result  from  a  Bayesian  approach. 
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IV.  IMPLEMENTATION  OF  THE  BAYESIAN  FORMULATION 


The  Bayesian  formulation  of  the  position  finding  problem  may  be  summarized  as 
follows.  The  relationships  which  characterize  the  site  geometry  f.  (x,  y)  and  the  joint 

density  of  the  observation  errors,  P(^).  are  assumed  known.  Prior  information  is  in¬ 
corporated  by  means  of  an  assigned  density  <y(x,  y).  The  principal  requirement  is  a 
feasible  computational  procedure  which  makes  use  of  these  known  or  assigned  quantities 
to  provide  an  estimate  of  the  position  of  the  source  of  radiation  corresponding  to  the  par¬ 
ticular  data  vector  yielded  by  the  observations.  This  estimate  is  to  be  chosen  so  that  the 
a  posteriori  risk  is  minimized. 


Bayes'  rule  provides  a  formal  statement  of  how  prior  knowledge  is  modified  by  ob¬ 
servations  so  any  algorithm  which  provides  the  a  posteriori  density  p(x,  y|  Q)  from  a 
knowledge  of  the  prior  density  er(x,  y)  and  the  data  vector  8  is  formally  equivalent  to  it. 
That  is,  the  algorithm  must  be  based  on  the  relationships  (3.2)  or  (3.2a),  namely: 


«)- 


y)  cr(x, 


3(91  *,  V)  cr(x,  y) 


x,  y)  tr(x,  y)  dx  dy 


We  note  that  the  conditional  density  p  [B  |x,  y)  is  readily  obtained  from  the  error  density 
p(\j/ )  which  is  assumed  to  be  known.  In  addition,  it  is  clear  that  if  the  density  function 
p(  8)  can  be  expressed  in  closed  form,  an  explicit  algorithm  for  effecting  the  mapping 
of  <r(x,  y)  into  p(x,  y|$)  is  assured.  The  derivation  of  an  explicit  algorithm  is  therefore 
seen  to  depend  on  obtaining  a  closed  form  expression  for 


x,  y)  o(x,  y)  dx  dy 

where  p(fflx,  y),  which  is  dependent  on  the  site  geometry,  has  a  complex  form  such  as 
that  given  for  a  specific  example  by  equation  (3.7). 

Rather  than  pursue  this  course  for  the  various  error  distributions  and  site  geometries 
that  could  exist  and  for  a  completely  arbitrary  choice  of  prior  density  we  adopt  an  ap¬ 
proach  based  on  the  following  argument.  Prior  knowledge  of  the  source  position  is  nec¬ 
essarily  imprecise.  Although  it  may  be  based  on  measurements  which  have  substantially 
greater  inherent  accuracy  than  those  under  discussion,  the  assumption  that  a  probability 
can  be  assigned  for  all  possible  values  of  the  coordinates  is  in  conflict  with  the  limited 
resolution  that  can  be  achieved  in  practice.  Consequently,  it  is  difficult  to  envision  a 
situation  where  all  of  the  prior  knowledge  cannot  be  accounted  for  by  assigning  proba¬ 
bilities  to  only  a  finite  set  of  source  coordinates.  This  leads  to  the  specification  of  a 
prior  density  having  the  form: 

0<x,  y)  =  ajk6(x  -  X.,  y  -  yk),  (4.1) 
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where  j  and  k  are  chosen  from  index  sets  of  finite  extent  such  as: 


j  =  -J.  -J  +  1 . -1,  0,  1 . J-l.  J 

k  =  -K,  -K  +  1 . -1,  0,  1 . K-l,  K 


and  6(x  -  x.,  y  -  y  )  denotes  a  generalized  bivariate  function,  usually  referred  to  as  a 

J 

delta  function,  whose  definition  and  properties  are  given  in  reference  (h).  Making  use  of 
these  properties  we  have, 


y)  <r(x,  y)  dx  dy 


=I  I  p(0|xr  yk)0jk 

j  k 


(4.2) 


which  is  in  closed  form.  The  a  posteriori  density  is  then  given  by  (3.2),  namely: 

p<*.  y|8>  -  El8|-’pfa)--ja-  (4-3> 


P (B  |x.,  yk)  oik  6(x  -  y.,  y  -  yk) 

21  P(0|x.,  yk)tTjk 

for  any  point  (x.,  y  ).  If,  for  particular  data  vector  $  which  results  from  the 

j  k 

observations,  we  write: 

Pjk  =  p(0l xr  yk> 


and  define 


qjk 


by 


the  relation: 


qjk 


-  Pjk  °jk 
2^pjk  °jk 


(4.4) 


(4.5) 


then,  on  substituting  in  (4.3)  we  obtain 

d(x,  y  1 8 )  =  qjk  6 (x  '  xj,  y  -  yk),  (4.6) 

for  each  value  of  j  and  k,  as  the  a  posteriori  density  corresponding  to  the  observed  data. 
Since  this  has  the  identical  form  to  that  given  by  (4. 1)  for  the  prior  density,  the  required 
algorithm  for  mapping  the  postulated  prior  density  into  its  a  posteriori  counterpart  is 
therefore  given  by  equation  (4.5)  above,  with  values  of  p.^  as  determined  by  equation (4. 4) 

from  a  knowledge  of  the  error  density,  site  geometry,  and  observed  data.  The  values  of 
p^  thus  reflect  the  specific  system  under  consideration  as  well  as  the  data  actually 
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recorded,  whereas  equation  (4.5)  is  a  restatement  of  Bayes'  Rule  applicable  to  the  form 
of  prior  density  chosen  to  describe  the  prior  knowledge. 

These  results  follow  from  the  assumption  that  all  of  the  prior  knowledge  may  be 
accounted  for  by  assigning  probabilities  to  a  finite  set  of  source  coordinates  and  do  not 
depend  on  how  this  set  is  chosen.  While  in  some  cases  a  rationale  for  selecting  this  set 
may  be  based  on  some  aspect  of  the  prior  knowledge,  it  will  often  suffice  to  allocate  one 
point  for  each  rectangular  area  enclosed  by  a  grid  system  spaced  in  accordance  with  a 
degree  of  resolution  that  is  considered  adequate  for  a  treatment  of  the  entire  problem. 
This  grid  system  assumes  a  comparable  role  in  the  formulation  of  the  a  posteriori 
density  giving  this  in  a  form  which  provides  a  convenient  starting  point  for  further  bear¬ 
ing  observations  or  for  implementing  discrete  search  procedures  such  as  those  described 
in  references  (i)  through  (k). 

Although  there  will  be  cases  where  the  overall  objective  is  merely  to  record  the 
estimated  source  position,  in  general,  some  subsequent  action,  such  as  search,  is  en¬ 
visioned.  This  action  should  be  based  on  the  a  posteriori  density  --  otherwise  all  of  the 
available  information  is  not  utilized.  Of  course,  once  thi6  density  is  known,  the  esti¬ 
mated  source  position  and  certain  related  information  is  readily  ascertained. 


For  a  quadratic  cost  function,  the  estimated  fix  is  obtained  from: 
x  =  E(x|$) 

=#  p  (x,  y  1 9)  dx  dy 

-II  V* 

j  k 

on  substituting  for  p(x,  yjg)  from  equation  (4.6),  and  from 


•II  A’jk- 

j  k 


(4.7a) 


Some  additional  properties  of  the  a  posteriori  density  which  may  prove  useful  can  be 
obtained  from  a  knowledge  of  the  moments  [1  ]  which  are  given  by: 

prs  =  E(xr  y8|0) 


-IP 


y  p(x,  y|0)  dx  dy 


II  *7  y?  ojk  • 

J  k 


(4.  8) 
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and  the  probability  that  the  source  is  located  in  any  given  region  may  be  determined  by 
summing  the  values  of  for  points  within  the  region. 

To  apply  these  results  to  the  example  given  in  Section  III  we  require  only  the  values 
p  appropriate  to  the  error  distribution  and  site  geometry  which  characterize  this  ex- 

ample.  These  are  readily  obtained  from  equation  (3.  7),  giving 


Pjk  =  p(0|  x.,  yR> 

n 

-i/§  —2 

/o  ■ 1  CP.  *  ; 


x  -n/2  /  x  i=l  i 

=  (2n)  (or..on)  e 


where 


x.  -  X. 

fi  (xj'  yk}  =  arctan  yt  -'V. 


and  the  observed  values  of  0  ,  0„ . 0  are  inserted.  The  appropriate  a  posteriori 

i  ^  n 

density  is  then  given  by  equation  (4.5)  and  the  estimated  position  of  the  source  by 
equations  (4.7)  and  (4.7a). 

We  observe  that  in  cases  such  as  this  where  the  observations  from  each  station  are 
independent,  it  is  not  necessary  to  await  the  reporting  of  all  n  bearings  before  commencing 
these  procedures.  A  simple  sequential  treatment  may  be  adopted.  Assume  the  rth 
station  to  be  the  first  to  report.  For  this  single  observation  the  appropriate  values  of  p  ^ 

are  given  by  the  univariate  conditional  density  p(©rjx,  y),  1  <  r  g  n.  Inserting  these 

values  in  equation  (4.5)  provides  the  a  posteriori  density  resulting  from  this  observation. 
When  the  second  observation  is  reported,  this  density  becomes  the  new  prior  density 
from  which  the  a  posteriori  density  corresponding  to  the  two  observations  is  computed 
using  the  same  procedure.  This  process  can  be  repeated  until  all  stations  have  reported 
or  until  sufficient  information  has  been  accumulated.  An  estimate  of  the  position  of  the 
source  is  readily  obtained  at  any  stage  by  means  of  equations  (4.7)  and  (4.7a). 

We  conclude  with  an  example  of  the  computation  for  the  special  case  referred  to  in 
Section  III  which  was  characterized  by  a  "flat-earth"  geometry,  normally  distributed 
errors,  and  independent  observations  from  each  station.  The  locations  of  the  stations, 
the  standard  deviations  of  the  error  distributions  and  the  recorded  observations  are 
taken  to  be: 


(miles) 


(miles) 


(degrees) 


(degrees) 

5.13 

36.7 

77.9 


It  is  assumed  that  the  source  of  the  radiation  is  shipborne  and  located  within  an  11  x  11 
square  mile  area  centered  about  the  origin  of  the  coordinate  system  and  that  this  area 
contains  a  small,  irregularly  shaped  land  mass.  In  addition,  the  position  of  the  source  is 
assumed  to  be  uniformly  distributed  over  the  residual  sea  area,  so  the  values  of  the  prior 
density  function  at  points  spaced  at  one  mile  intervals  are  as  shown  in  table  I  from  which 
the  location  and  shape  of  the  land  mass  is  readily  apparent. 

Since  the  procedure  for  computing  the  required  a  posteriori  density  is  repeated  for 
each  point  of  the  11x11  array  resulting  from  this  choice  of  spacing,  we  consider  an 
arbitrary  point  (x.,  y, )  at  which  the  value  of  the  prior  density  is  a.,  6(x  -  x.,  y  -  y, ). 

J  K  ]K  J  K 

For  this  point,  the  true  values  of  the  bearings  from  each  station  are  given  by  equation 
(3.4),  namely 

x.  -  X. 

a.  =  f.  (x.,  y  )  =  arctan  — 1 ^ - 

1  1  1  V  yk  -  Y. 

for  i  =  1,  2,  3.  The  numerical  value  of 

Pjk  =  P(8]xr  yk) 

is  then  obtained  from  equation  (3.7).  This  is  repeated  for  each  point  to  give  the  array 
shown  in  table  II.  From  this,  the  required  a  posteriori  density  is  computed  by  repeated 
application  of  equation  (4.5)  and  the  estimated  source  position  by  means  of  equations  (4.7) 
and  (4.7a).  The  numerical  values  so  obtained  are  given  in  table  III. 

This  example  is  reconsidered  for  the  case  in  which  the  three  observations  are  not 
simultaneously  available  by  means  of  the  sequential  procedure  outlined  above.  Table  IV 
gives  the  values  of  the  density  function  prior  to  any  station  reporting  and  table  V  the 
a  Dosteriori  density  based  on  the  single  observation  reported  by  Station  1.  This  density 
assumes  the  role  of  a  new  prior  density  function  for  the  computation  of  the  a  posteriori 
density  corresponding  to  the  second  observation  to  be  reported.  Table  VI  gives  the 
numerical  values  at  this  stage  and  a  similar  computation  for  the  observation  reported  by 
Station  3  yields  the  final  a  posteriori  density  corresponding  to  all  of  the  observations 
shown  in  table  VII.  As  would  be  expected,  this  is  identical  to  table  III.  The  estimated 
position  of  the  source  at  each  stage  is  readily  computed  from  the  corresponding 
a  posteriori  density. 
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TABLE  I 


PRIOR  DENSITY  FUNCTION 


o :  j  s 

0.00673 

C  .  00803 

0 . 0  0  0  <1 3 

'*0.000  3  3 

0.00803 

0.00893 

0 . 0  0  p  9  9 

0  ,  c  0  P  9 1 

0 .00893 

0 .008  *3 

"  "■ a  *  i 

0.009  >  3 

.  008  *3 

0. 008 >3 

0  .  O  0  9 1 3 

0 . 0  O  R  9  3 

0.00893 

0.09893 

0 . 0  C  H  9  3 

0.00893 

0.00893 

0 . 0  0  P  3  3 

.  C 080  3 

0.00803 

0.00893 

0.00893 

0  .  OOm93 

9.00  P  '>  3 

0.00891 

0 .00893 

0.00893 

0  0  5  >3 

0 . 008  S3 

0.00833 

0.00313 

0.00833 

0.00893 

0.00893 

0 . OOP  9  3 

0.00893 

0.00893 

0. 008°3 

D  33  13 

0 . 00811 

3.00803 

0. 00803 

0. 008-13 

0.00893 

0.00893 

0 .90393 

9.00893 

0.00893 

0  .  00893 

0  38  t  3 

0.0C313 

0.00333 

0 . 008  j 3 

0.00813 

0.908  93 

C  .008-3 

0.00? 33 

9 . 0  0  P,  1  3 

0 . 00B93 

0 .00893 

00313 

0 . 0C8  33 

0.00893 

0. 0  03 --3 

o.oo.-;  9  3 

0.00093 

C  .00893 

0 

0 

0 .00093 

0 .00993 

0  C  8  1  3 

0 . 003 >3 

0.03333 

0.008:3 

0.00903 

0.00893 

0 

n 

0 

0 .00893 

0 .00893 

0  0933 

0.00033 

0.00893 

0.00303 

0.00393 

0.00803 

0 

0 

C 

0.00893 

0.00893 

0  0  3  *  3 

o .  o  c  a  i  3 

0.00603 

0.00813 

C  .  008  3  3 

0.00893 

0 . 00099 

0 

0.00893 

C  .  00893 

0.00893 

3  05  >3 

0.00133 

0.008-0  3 

0.00833 

0. 00a»3 

0. 0  08') 3 

0.00B93 

0.  101  »3 

0.00893 

0.00893 

0.00893 

7.91  r  *■> 

F.  .53 r’s 

1.28  .-■*• 

2 . 1  l.p”‘l 

3 . 9  7F~i 

8  .  5 

1.65p~4 

2.03 F~u 

5 .  ?if”r 

1  .  1  3  F~ l* 

1 

4. 

r~ 

0 

3 . 1 6F”'» 

••  .  -pf-* 

1  .u3  r-*. 

2.5ir  m 

3 . 6BF% 

-  .  i  3  F "  \ 

1  . 73F~4 

2 . 3  3  r**u 

u. 13F”u 

: . 

2 .0“  r-** 

3 .  air*1* 

4  47£'”u 

<  •*  7  r*".. 

: .  3ir“i« 

3  .f  ?.r% 

4.3  6.p*4 

1.3  - 

2. so "'4 

3 . 77r"w 

4 . 0  F  4 

1  .  r.  IF”- 

5 .  ■■  5  r“- 

3 . 7 IF* 4 

4.47F"- 

1  .  c  •  “  ”  - 

2 . 6  7  5  ”  4 

3  .  C  5F~4 

4 . 1  3  F  ”  4 

1  .6  r  "_4 

2 . 5  P  "”*• 

3 , 3  7  F  ”  4 

3.«-  IF” 4 

TABLE  II 

VALUES  OF  p.k 


2  .  *lci.p*4 

3. 5  IF” 4 

3.r.sr*4 

3  , 

3 ,53p“4 

•* .  01F_4 

3 . 8  7F~4 

3, 

4 . 07.p”*. 

4 . 4  3F~4 

u  .  0  8F~4 

3  , 

4 . 5  4. <*“4 

4 . 7lF_4 

4.13F"4 

1 

4  .  C  7  F“  4 

4.8?F“4 

u  .  32F”‘v 

7 

5 . 94/7 ~4 

4 . 7  4  F  ”  4 

3.75r"4 

2 

5 . 9  0  :’~4 

4 .43F"4 

1 . 35F"4 

0 

4 . 7.'.r~4 

4 . Q  4  F  ”  4 

1  ,  P7F~4 

1 

4 . 3  4  r — 

3 . 49F~4 

7 . 34F_4 

1 

3 . 7  9  r”  4 

?.98F-4 

1 .92F“4 

3 . :  4  f”  4 

?  .  ? OF "4 

1 . 34.r'"4 

f 

07F“4 

2 .27F_4 

1  . 44F”4 

7.9 2F*  5 

1  OF*  4 

2 .24F*4 

1 .35F"4 

7  . OOF'S 

1  OF ‘4 

0 . 1 2F~4 

1  .21F‘4 

5 . 9  2  f” 5 

Of  F“‘» 

1  .  9  3  F~  4 

1  .04  F~ 4 

4.8  0F~  5 

,0  3  .**”'» 

1  .f  °r~4 

e.56F”5 

3 . 7 1F~  5 

.  5CF"4 

1 .41F"4 

r . 7 1 r~  5 

2 . 73F”5 

.  1  1F*4 

*  .l?.p“4 

5 .  0  2  F  ”  6 

1 . 9  2  F” 5 

.79F"4 

R  .i»8F"5 

3 . 5  7  F” 5 

1 . 27  F” 5 

.  3  1 F  ~  4 

9.1  IF"  5 

2.“ IF” 5 

8 . 02r”e 

.  r>  3F**5 

4.1 SF”5 

1 . 5UF"5 

4 . 78F“6 

.  (  n-“~r. 

?  .  7  1 F  "  5 

9 .29"”' 

2 .  c  sf”c 

TABLE  III 

A  POSTERIORI  DENSITY  FUNCTION 


0 . U  *  1 05 

0.00245 

0 . 0  C  4  7  9 

0 . 00791 

0.9111 

0.0131 

0.0133 

0.0115 

0.00051 

0.00541 

0 .00297 

0.00146 

0 . 0  r  3  2  P 

0.00613 

0.00384 

0.0132 

0.015 

0 .0  145 

0.0119 

0.00039 

0.00506 

0 . 00262 

0.00195 

0.00425 

0 . 00775 

0 . 7 1 1  fl 

0. 01  S3 

0.0166 

0.0153 

0 . 0  1  1  9 

0 .00796 

0 . 00454 

0.00222 

0.00254 

0.00534 

0  .  00938 

0.9138 

0.017 

0.0170 

0.0155 

0.0115 

0.09724 

C.0039 

0.0018 

9.0032 

9 . 0  0t  5 

0.011 

0.0155 

u  .0183 

0.0181 

0.015 

0.0106 

C  .  0003 1 

0 .003  2 

0. 00139 

0 . 90387 

C  .  0  0  7  Fj  3 

0.0124 

0.0168 

0.016  » 

O.C 1 78 

0.014 

0.00935 

0.09526 

0 .00251 

0 .00102 

~ . C  94  5  7 

0.00865 

0.0135 

0.0174 

0. 9187 

0.0168 

0.0126 

0 

0 

0.00188 

0 . 000717 

0 . 09513 

0 . 009**  1* 

0  .0141 

C  ,  0174 

0.017? 

0.0151 

0 

0 

0 

0.00134 

0 . 000477 

9.0956ft 

C  .  00  *7  1 

0  .  C  1  4  7 

0.0167 

0.01*  3 

0.0131 

0 

0 

c 

0 .000902 

0.000301 

9.00594 

0  .  C0993 

C.C137 

9.0154 

9.0142 

0. 010P 

0.0068 

0 

0.00157 

0.000576 

0 . 00017° 

5.00599 

9 .0096  5 

0 . 7  1  2  6 

0.0135 

0.0118 

n . no84f 

0.00602 

0.00247 

9 . 0  0  1  0 1 

0.900348 

0 .0001 

x  =  -  0.78 
y  =  0.24 
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TABLE  IV 

PRIOR  DENSITY  FUNCTION 


r .  o  o  B  ■»  a 

0.00803 

0.00803 

0.00893 

0. 00R03 

0.  00893 

0 . 00P93 

0 . 0  08*11 

0.00891 

0.90 f 9 9 

'  a 

r- .  o  o  p  9  a 

0.00803 

0.00803 

0, O0893 

0.00893 

0.00893 

0.00893 

0.00893 

0.90891 

0.00  <*  9  9 

f  '  - 

P. 00093 

0 . 0  OP  0  3 

0.00803 

0.00893 

0.00893 

0.00893 

0. 00893 

0 . 00  8  9  3 

0.00891 

9 . 0  0  8  "  1 

r  .  0  9  ?■ 

91 

0.00833 

0.00003 

0.00803 

0.00803 

0.00893 

0.00693 

0.00803 

C  .  0  0  F.  9 1 

0 . 0.08*1 3 

0 . 9  0  Y.  1 9 

1. OOR93 

0.00803 

0.00893 

0.00803 

0.00893 

0.00093 

0.00093 

0 . 00R9? 

0 . 0  0893 

0.9089  9 

0.90  •: 

,■5 

0.008  13 

0.00803 

0.00803 

0.00803 

0.00893 

0.00893 

0.00893 

0.00891 

0.00893 

0.00891 

G.  QC'i 

9  9 

0 . OOP  0  3 

0. 00803 

0.00803 

0.00893 

0.00893 

0.00C93 

0.00893 

0 

0 

0.00893 

0 . 0  0  8 

'■3 

0 . 00803 

0.00803 

0.00803 

0.00803 

0.00813 

0.00893 

0 

0 

0 

Ci .  0  C  ft  3  9 

0 . 0  06 

.  3 

0 . 00803 

0.00803 

0.00803 

0. 00893 

0.00893 

0.00893 

0 

0 

9 

0  .  OOP 93 

.  9  0 

0.00803 

0.00803 

0.00803 

0.00893 

0 .00893 

0.00893 

0.00893 

0 

0,00891 

Ci  .  08  •  9 

j  3 

0.00833 

0.00883 

0.00803 

0 . 00893 

0.00893 

0. 00893 

0.00893 

o .  o  o  o  i 

0.90899 

0 . 0  08  9  • 

0.0108 

0.0137 

0 .0156 

TABLE  V 

A  POSTERIORI  DENSITY  FUNCTION  (1) 

0.0158  0.0143  C.0117  0.00848  0.f''5r.4 

'.''332*. 

c .  *  o  i  ? : 

0.011 

0.0139 

0.0157 

0.0157 

0.0141 

0.0113 

0.00807 

0.19537 

0  0*1  2 01 

C. 0112 

0.0141 

0.0157 

0.0156 

0.0330 

0.0109 

0.00766 

6 .004* 

n  <■>  n  ?  r 

0.0142 

0.0158 

0.03  55 

0.0136 

0.0105 

0.00724 

0.0144 

0.0158 

0.0154 

0.0133 

0.0101 

0 .  °  0  r,  ft  ? 

C.O04  0'* 

o  .  f  c ;  ' 

0.0118 

0.0146 

0.0159 

0.0153 

0.0129 

o. ooor c 

0.0064 

0. 00275 

0.r 0195 

0  .  r.  0  0  :■ 

i-  r  f, 

0.01? 

0.0147 

0.0153 

0.0151 

0.0126 

0.00925 

0 . 0  0  r.  9  1 

n 

o 

0.012? 

0.014? 

0.0150 

0. 0149 

0.0122 

o.  coen 

0 

9 

0 

0.0124 

0.015 

0.0159 

0.0147 

0.0119 

0.00826 

r, 

0 

0 

'  f-r,f  -  r  ■) 

C . 01 26 

0.0152 

0.0159 

0.0145 

0.0115 

0.00731 

0 . 00477 

0 

o .  r  c  1 1 1 

0 . 0  C  0  4  7  1 

0.0153 

0.0158 

0.0142 

0.011 

0.00745 

0.00430 

0.00224 

0.PG1 

0.000352 

0.00  0 1  1 1 

x  =  -2.1 

y  =  -0.25 

0.000566  0.00132 

0.0026 

TABLE  VI 

A  POSTERIORI  DENSITY  FUNCTION  (2) 

0.0043?  0.00603  0.00728  0.00743  0.00647 

0.00483 

0.0031 

0 . 0  C  3  7  2 

0. OOOBOS 

C. CO104 

0.0035 

0.00562 

0.00762 

0.00876 

0.00855 

0.00731 

0.00506 

0.00308 

0.  CC1G  2 

0.00115 

0.00253 

0.00466 

0.00722 

0.00942 

0.0104 

0.00967 

0.00766 

0.0C518 

0.00299 

0.00148 

0.00161 

0.00345 

0.00614 

0.00916 

0.0115 

0.0121 

0.0107 

0.00809 

0.00518 

0.00283 

0. 00133 

0.00225 

0.00465 

0.00798 

0 . 0  1 1 4 

0.0137 

0.0138 

0.0117 

0.00835 

0.00507 

0.00261 

0 .00115 

0.0031 

0.00619 

0.0102 

0.0141 

0.0161 

0.0155 

0.0124 

0.00842 

0.004R2 

0.00234 

0  .  UCC97 

0.00422 

0.00812 

0.0129 

0.017 

0.0186 

0.017 

0.0129 

0 

0 

0.00204 

0 . 0  **0  7  9 

0 . 00c67 

0 .0105 

0.016 

0.0202 

0.021 

0.0182 

0 

0 

0 

0.00172 

0. 0006? 

0.00751 

0.0134 

0.0195 

0.0235 

0.0232 

0.019 

0 

0 

0 

0.0014 

0.00047 

0.00379 

0.0167 

0.0234 

0.0267 

0.025 

0.0194 

0.0124 

0 

0,00295 

o  .non 

0  .  0003* 

0.0125 

0.0205 

0.0273 

0.0296 

0.0263 

0.0192 

0.011C 

0.00577 

0 .0024 

0.000836 

0.00024 

x  =  -1. 

y  =  -i- 

0.00105 

1 

1 

0 .00245 

0.00473 

TABLE  VII 

A  POSTERIORI  DENSITY  FUNCTION  (3) 

0.30791  0.C111  0.0131  0.0333  0.0115 

0  .  P  0  8  f>  1 

O.PPf.u 1 

0.00146 

0 . 00328 

0.00619 

0.009P4 

0.013? 

0.015 

0.0145 

0.0119 

0  .  P  C  8  1 9 

9  .  ?  C  5  P  r 

■  v  . 

0.00195 

0.00425 

0.00775 

O.Ollfl 

0.0153 

0.0166 

0.0153 

C.01  19 

0  .  P  n  7  n  f. 

n  .  p  p  4  *  4 

,  r  r  7  " 

0.00254 

0.00534 

0.00938 

0.0130 

0.017 

0.0176 

0.0155 

0.0115 

0.00724 

0 . 0  a  1  ;< 

o .  c :  i  >• 

0.0032 

0.0065 

0.011 

0.0155 

0.0103 

0. 01R1 

0.015 

p.mor 

0  .  p  0  0  3  1 

P  .  "  c.  ?  .» 

0.00389 

0.00763 

0. 0124 

C.O1G0 

C . 01 09 

0.0178 

0.014 

0 . 0  0  r'  3  s 

n  .  1 0  5  2  f 

C  1 

'  .  '  ~  1  v . 

0 .00457 

0.00865 

0.0135 

0.0174 

0.0187 

0.0160 

0. 01 20 

0 

0 

0  .  P  r 

0 .00519 

0.00944 

0.0141 

0.0174 

0.0178 

0.0151 

0 

0 

0 

or  l  in 

o  u  - 

0 . 00566 

0.00091 

0.0142 

0.0167 

0.0163 

0.0131 

c 

0 

0 

n .  n  t)  n  i  o  2 

C  .  •;  1 i 

0.00594 

0.00999 

0.0137 

0.0154 

0.0142 

0.0108 

0.00(1! 

0 

0.10157 

r. .  i  ?  r  *  7 1 

2  .06017 

0.00599 

0.00965 

0.0126 

0.0135 

0.0118 

0.00846 

0 .005  0? 

0 . 0  0  2  4  7 

P  .  0  1 1  0  1 

P  .  1  0  •  •  9  •*  P 

•'  .  0  0  f]  i 

x  =  -0.78 
y  =  0.24 
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V.  SUMMARY  AND  CONCLUSIONS 

In  this  communication  the  position  finding  problem  is  considered  from  a  standpoint 
which  includes  the  following  features: 

(a)  The  position  of  the  source  of  radiation  is  assumed  to  act  as  a  random  variable 
for  which  a  prior  density  can  be  assigned. 

(b)  The  position  finding  system  is  assumed  to  consist  of  a  set  of  geometrical 
relationships  from  which,  in  the  absence  of  uncertainty,  the  position  of  the  source  could 
be  determined,  together  with  an  error  distribution  to  account  for  the  uncertainty. 

(c)  Emphasis  is  placed  on  obtaining  closed  form  expressions  for  the  decision  pro¬ 
cedures  established  as  a  means  of  avoiding  the  uncertain  accuracy  of  approximate  pro¬ 
cedures  and  the  computational  complexity  of  iterative  procedures. 

The  problem  is  recognized  to  be  fundamentally  one  of  ascertaining  the  manner  in 
which  prior  knowledge  is  modified  by  observed  data  and  the  proposed  solution  is  based 
on  a  Bayesian  treatment  which  assumes  a  particular  formulation  of  the  prior  knowledge. 
The  formulation  is  well  suited  to  the  frequently  occurring  case  where  this  knowledge  is 
empirical  in  nature  and  is  nonrestrictive  in  the  sense  that  any  prior  density  can  be 
approximated  with  arbitrary  precision.  This  approach  yields  a  simple  algorithm  for 
providing  the  required  a  posteriori  knowledge  in  a  similar  formulation.  This,  in  turn, 
characterizes  the  uncertainty  in  the  presence  of  which  any  subsequent  action  is  initiated 
in  a  form  that  is  particularly  convenient  where  this  action  consists  of  making  further 
bearing  observations  or  implementing  certain  discrete  search  procedures.  The  required 
estimator  of  the  source  position  is  also  readily  ascertained  by  means  of  this  approach. 

When  compared  with  procedures  in  current  use  the  Bayesian  treatment  proposed  here 
offers  the  following  advantages: 

(a)  It  provides  for  the  incorporation  of  prior  knowledge. 

(b)  It  provides,  in  addition  to  an  estimate  of  the  source  position,  all  relevant 
a  posteriori  knowledge. 

(c)  It  provides  simple  computational  procedures  for  obtaining  these  results. 

(d)  It  is  applicable  to  a  diverse  range  of  position  finding  systems. 
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